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Abstract The joint use of high-resolution data from PROBA2 and SDO satel- 
lites and LASCO/SOHO coronographs enabled us to examine early stages of 
initiation and propagation of six limb CMEs registered in June 2010 - June 2011. 
For five events under consideration, the CME initiation is marked by filament 
(prominence) eruption or by a loop-like structure having another nature. Subse- 
quently, several loop-like structures having higher brightness and following each 
other at different velocities appear in the region of the CME initiation. The CME 
frontal structure is formed by these loop-like structures. The time-dependent 
velocities and acceleration of the ejection front have been obtained for all CMEs 
under consideration. We have drawn a conclusion about the possible existence 
of two CME types dependent on the time profile of their velocity. The first 
CME type comprises the ejections whose velocity decreases abruptly by more 
than 100 km/s after having reached the maximum; it thereupon passes to slow 
deceleration. The second CME type is formed by the ejections whose velocity 
varies insignificantly after reaching the maximum. The CME angular size is 
shown to increase up to threefold at the initial stage of propagation; it increases 
twofold 3.5-11 minutes after the first measurement of this parameter. When 
considering 3 CMEs, we see that their broadening exceeds their extension in the 
longitudinal direction during a certain period of time at the initial propagation 
stage. 

Keywords: Coronal Mass Ejections, Initiation and Propagation; Corona, Struc- 
tures 



1. Introduction 



Physical mechanisms of generation of coronal mass ejections (CMEs) still re- 



main unclear in many respects (Forbes et al., 2006). This burning issue can 



be addressed due to the use of data with high temporal and spatial resolution 
that will allow us to study initiation and initial propagation stage of CMEs. In 



2010, Solar Dynamics Observatory (SDO) (Pesnell, Thompson, and Chamberlin, 
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2012) was launched with the on-board Atmospheric Imaging Assembly (AIA) 
instrument (Lemen, Title, and Akin, 2012) having unique characteristics. AIA 



enables us to observe the Sun in the wavelength range from extreme ultraviolet 
(A = 9.4 nm ) to ultraviolet (A = 170.0 nm ) with high temporal (10 seconds) 
and spatial resolutions (0.6 arcsec pixels). Consequently, it will be possible to 
reveal regularities of the CME generation and the initial phase of the CME 
propagation. 

The additional use of data from the PROBA2 (Project for On-Board Au- 
tonomy) satellite launched in 2009 will contribute to this study, too (http: 
//proba2. oma.be/index. html/about/). The on-board instrument SWAP 
(Katsiyanni, Berghman, and Swap Consortium, 2010) observes the Sun in ex- 



treme ultraviolet at 17.4 nm, whereas LYRA (http://proba2.oma.be/index. 
science/) monitors the solar irradiance in the wavelength range from soft X-rays 
to ultraviolet. Though the temporal (1 minute) and spatial (3.17") resolutions 
of SWAP are worse than those of AIA instruments, it has a big advantage: its 
field of view is significantly larger (54 arcmin as against 41 arcmin (SDO)). 

Establishment of experimental regularities of the CME propagation imme- 
diately after its initiation is of importance when studying CME features. Over 
the last decade, several studies were devoted to different aspects of the initial 
stage of the CME propagation, using data from TRACE, SOHO, STEREO, 



GOES satellites (see papers ( Gallagher, Lawrence, and Dennis, 2003 Zhang and 



Dere, 2006 Maricic, Vrsnak, and Rosa, 2009[ )). First of all, kinematics of the 
CME propagation was examined. Three periods of the initial propagation of 



CMEs were established (Zhang and Dere, 2006): slow increase in velocity within 



several tens of minutes, rapid increase in velocity (the main CME acceleration) 
within the period from several minutes to several tens of minutes (sometimes up 
to several hours) when the acceleration amplitude reaches the value of several 
km/ s 2 in some ejections. The third period is the propagation phase of CMEs (i.e., 
the 'quiet' propagation of CMEs with an acceleration that is small or moderate 
in absolute magnitude). 

A close connection of the CME acceleration with energy release was estab- 



fished in the CME-related solar flares ( 


Zhang and Dere, 2006 


Maricic et al, 


2007 


Temmer et al, 2008 


Temmer et al, 2010 


). There was a positive correla- 



tion between the duration of the main CME acceleration tACC and the time 
of increase tSXR in intensity of soft X-rays ISXR from the CME-related flare 
area up to maximum values ( Zhang and Dere, 2006| Maricic et al, 2007). The 
CME acceleration a(t) and hard X-ray emission Ijix(t) ( or its related derivative 



dlsxFt/dt) from the flare area were shown to be often synchronised (Maricic 
et al., 2007} Temmer et al., 2008 |Temmer et al, 2010[ ). An inverse correlation 



dependence of acceleration V m ax/t a cc on t acc (or V ma x/tsxR and tsxn) was 



established ( Zhang and Dere, 20061 Vrsnak et al., 20071. 



Due to unique capabilities of the instruments on board SDO, new results 
on the CME initiation and its initial propagation stage have been obtained 



(Patsourakos, Vourlidas, and Stenborg, 2010). With a CME as an example, we 



have shown that the CME initiation has three stages: the stage of slow self- 
similar expansion, short-term strong hyperexpansion in the transverse direction, 
and one more stage of self-similar expansion in the SDO field of view. There is 
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also a strong short-term CME acceleration with a subsequent deceleration. Note 
that the conclusion about the existence of a period of the CME hyperexpansion 



in the transverse direction was drawn from PROBA2 data later in (Fainshtein 



and Egorov, 2011 ). 



In this paper, we go on studying the CME initiation during the new solar cycle, 
using SDO/AIA data. We also examine the initial stage of the CME propagation 
with the use of data from the SDO, PROBA2 and SOHO spacecraft. The goal of 
this study is to get new experimental results concerning the CME kinematics and 
time variations of their geometric parameters that could be used to construct 
adequate models of the CME generation and propagation at the initial stage. 



2. Data and investigation techniques 

We have analysed CMEs of the new solar cycle, recorded during the period 
from June 2010 to June 2011. The events to analyse had to meet the following 
requirements: 

1. two spacecraft - SDO and PROBA2 - must observe CMEs. Resolution of the 
corresponding CME images to analyse must be sufficiently high. In particular, 
at least 4 consecutive images must represent each event in order to enable 
position of the CME front be determined and time-dependent velocity and 
CME acceleration to be calculated; 

2. events must be related to X-ray flares; 

3. events must be limb or close to limb ones; 

4. each ejection selected using SDO and PROBA2 data must correspond to a 
CME recorded in the LASCO field of view, according to the catalogue data . 

There were 9 events that met these requirements. CME observed on 18 August 
2010 was not considered, since its initiation and propagation were complex. We 
did not consider the non-limb event of 14 August 2010 and the event of 24 
February 2011 either, since interpretation of its initial data still goes on. An event 
was considered to be limb, if the CME-related flare was at a distance > 60° from 
the centre of the solar disc (in the heliographic coordinate system) . There were 
5 such events: on 13 May 2010, 11 February 2011, 08 March 2011, 27 April 2011, 
and 07 June 2011. The flare position in 2011 was determined as the position of 
the source of X-ray emission in the range from 6 to 12 keV, according to Reuven 
Ramaty High-Energy Solar Spectroscopic Imager (RHESSI) data (http://sprg. 
ssl.berkeley.edu/~tohban/browser/). The event of 27 March 2011 was difficult to 
relate to a flare. Its source was probably just behind the limb. We estimated 
latitude of this CME's source: it was « 25° — 35° relative to the limb. 

The CME initiation and propagation at the initial stage have been studied 
using SDO/AIA data and FITS images of the Sun (http:/ /jsoc. stanford.edu/ajax/ 
lookdata.html), data processing level 1 in EUV lines A = 17.1 nm and 21.1 nm 
that arc formed when plasma temperature is ~ 9 105 K and 2 106 K. To analyse 
the CME kinematics, we used calibrated solar images (http://proba2.oma.be/ 
swap/data/bsd) taken by SWAP in the emission line having a wavelength of 17.4 
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nm. To determine CMEs using SDO/AIA data, difference images were plotted 
for each pixel with the use of the following relation: 

A7=(7 2 ^-7 1 ^)/7 

exp2 -L expl 

Here, I2 is the image brightness at a later moment of time; 7i, at one of the 
previous moments of time; 7o, at the third moment of time which is earlier than 
Ii . We took account of the fact that exposure time T exp i , T exp2 and T exp0 can be 
different when obtaining three images. To determine CMEs from SWAP data, 
we used only difference A7 = I2 — Ii ■ 

Data (http://sharpp.nrl.navy.mil/cgi-bin/swdbi/lasco/images/form) from Large 
Angle and Spectrometric Coronagraphs (LASCO) aboard SOHO were used to 
examine the CME propagation at distances R > 2i?o (R is counted from the 
centre of the solar disc in the plane of the sky, _? is the solar radius). To 
get information about soft X-ray emission, we used data from Geostationary 
Operational Environmental Satellites (GOES) (http://goes.ngdc.noaa.gov/data) 
and PR0BA2/LYRA (http://proba2.oma.be/index.html/sciencej). For one of the 
events (on 08 March 2011), we compared the time variation in the main accel- 



eration with the time dependence of hard X-rays recorded by RHESSI (http 
//sprg. ssl.berkeley.edu/~tohban/browser/) over the energy range from 25 to 50 
KeV. 

Information about eruptive filaments (prominence) from the closest proximity 
to the CME generation region was provided by the telescopes observing the Sun 
in the H-alpha line (http://swrl.njit.edu/ghn_web/) and by Coronado SolarMax 



H-alpha Imager at Mauna Loa Solar Observatory (MLSO) ( http:/ /mlso.hao.ucar. 
|edu/mlso-dataXORONADO-2010.htmi] ). 

The CME velocity Vi was calculated using formula 

Vi = Lr+1 Lf 
ti+i — u 

where Lf is the distance from the CME front to the chosen point on the visible 
solar disc, generally along the curved trajectory, at moment of time ti. We sup- 
posed that velocity V% was at the moment tf +l — + U)/2. Obtained values 
Vi(t) were approximated by B-splines, taking account of the error in determining 
Vi(t) at each moment of time. This resulted in the smooth function V[t). In the 
field of view of LASCO coronagraphs, we determined the HCME front position 
from the catalogue (http://cdaw.gsfc.nasa.gov/CME_list). The CME acceleration 
was determined from ratio a(t) — dV{t) /dt. The minimum error in calculating 
velocity was determined from ratio AV m i n = 25P '/ '(ti+i — ti), where 5P is the 
spatial resolution of the instrument. The CME front position was determined 
visually from solar images, whereas the error in calculating velocity was deter- 
mined as root-mean-square deviation of the calculated CME-front velocity from 
the mean value during N measurements (N = 5 - 10). We supposed that the 
CME front position could be determined accurate within 0.5i? in the field of 



view of LASCO at R > 5Rq (Shanmugaraju et al, 20101 
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To determine CME velocity V(t), we used not only values of the CME front 
velocity obtained using data from a certain instrument, but also 'intermediate' 
points V(T 1>2 ) = {R{T 2 ) - R(Ti))/(t T2 ~ t Tl ), where R(T 2 ) is the CME front 
position at moment tr 2 of its first registration by telescope 2, R(Ti) is the CME 
front position at moment tr 1 of its last registration by telescope 1. For instance, 
T x - SDO/AIA, T 2 - LASCO C2. 

To determine reliably the eruptive prominence (flux ropes, according to our 
terminology), loop-like structures that form CMEs, and the CME frontal struc- 
ture, special films were made, which allowed us to determine the structures that 
moved steadily. In what follows, we will use the term 'flux rope' to denote loop- 
like structures from hot coronal matter. To determine the magnetoplasma flux 
rope that preceded the eruption onset several tens of minutes and that would 
later form a CME, we used source images of the Sun after special processing. 
Numerous coronal loops were observed at these images. The relevant loop was at 
the place where it would be later observed as an eruptive flux rope on difference 
images. 

3. Results 

Fig. 1 exemplifies CMEs on 07 June 2011, according to observations of solar 
radiation in wavelengths of 17.1 and 21.1 nm (SDO/AIA data), and CMEs 
recorded by SWAP/PROBA2 and LASCO C2 coronagraph. Analysis of all solar 
images available from SDO/AIA instruments has revealed that initiation of this 
CME is due to motion of a small structure (see arrows 1 in Fig. 1). As we shall 
see later, this structure is an eruptive filament (prominence). When CME is 
under formation, the eruptive filament is preceded by numerous motionless and 
moving loop-like structures that show an increased brightness. In two lines, sets 
of such loops, however, are observed in a different way, and they vary over time 
in a different way. The moving loops have different velocities, and, as we shall 
see subsequently, some inner loops move at higher velocities than the outer ones. 
The inner loops finally catch up with the outer ones, thus leading to formation 
of a frontal structure (i.e., the outer region of the formed CME, showing an 
increased brightness). Let us compare the CME frontal structure in the two 
lines at the moment of time t =06:26:48. In the line at 17.1 nm, the frontal 
structure is closer to the ejection source than that in the line at 21.1 nm. The 
brightest part of the frontal structure is in the line at 17.1 nm (see Fig.l (c4)). 

CME of 07 June 2011 is related to the filament eruption, which is con- 
firmed by Fig. 2. It illustrates the filament (prominence) eruption from the 
region of the CME initiation, according to H-alpha line observations (data from 
Kanzelhohe Solar Observatory http://www.kso.ac.at/index_en.php) (there are no 
H-alpha data during the onset of the filament eruption). Fig. 2(d) shows the 
eruptive prominence observed in the H-alpha line and superposed onto the loop 
recorded in the line at 17.1 nm. Similarity in location of structures in different 
lines confirms that the moving line observed in extreme ultraviolet is the eruptive 
filament. 

Fig. 3(a) presents time dependence of velocity of the eruptive filament (EF) 
and of two loop-like CME structures that occurred after initiation of the EF 
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Figure 1. The event of 07 June 2011. al - a5 are the difference solar images obtained using 
SDO/AIA data in the spectral line A = 17.1 nm; bl - b5, in the lines A = 21.1 nm; cl - c2 are 
the difference solar images obtained using PROBA2/SWAP data in the line A = 17.4 nm; c3 
is the difference image of the solar corona, according to LASCO C2 data; c4 are the brightness 
scans along segments of the straight black line presented in Fig. a5 and b5. Solid line in Fig. 
c4 denotes the brightness scan in Fig. a5; dashed line, in Fig. a6. Arrows No. 1 denote the 
eruptive prominence; arrows No. 2, the inner loop when CME is under formation; arrow No. 
3, the outer loop; arrow No. 4, the formed frontal structure. Each image is provided with a 
coordinate axis allowing us to determine location of any element on the image, relative to the 
centre of the solar disc, in terms of seconds of arc. In Fig. a5 and b5, dashes perpendicular 
to segments, along which brightness was scanned, denote inner and outer boundaries of the 
region of the most significant variation in brightness in the CME frontal structure. In Fig. c4, 
they are marked off by vertical lines: solid ones for the image in the line A = 17.4 nm; dashed 
ones, in the line A = 21.1 nm. 
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Figure 2. (al - a3) is the eruption of filament related to CME of 07 June 2011, according 
to measurements in the H-alpha line. The white line indicates the boundary of the eruptive 
filament. a4 is the comparison between eruptive filaments in the H-alpha line (white dashed 
line) and in the line A = 17.1 nm (white solid line) at 06:37:02 (for the line at 17.1 nm) and 
at 06:37:04 (for the H-alpha line). Coordinate values on the coordinate axes are in terms of 
seconds of arc from the centre of the solar disc. 
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Figure 3. (a) is the time dependence of velocity of the eruptive prominence (solid line) and 
of two loop-like CME structures (dotted line denotes the outer structure; dash-and-dot line, 
the inner one), according to SDO/AIA data for the event of 07 June 2011; (b) is the time 
dependence of acceleration of the said structures; (c) is velocity V(t) obtained from joint 
SDO/AIA (squares), PROBA2/SWAP (triangles) and SOHO/LASCO (circles) data - solid 
line, and time-dependent soft X-ray intensity ISXR(t) - dashed line; (d) is the acceleration of 
the CME frontal structure a(t) - solid line - and time dependence dISXR/dt - dashed line. 



motion. Fig. 3(b) demonstrates time profiles of their acceleration. These loop- 
like structures (inner and outer ones) are marked off by arrows 2 and 3 in Fig. 
f. Results obtained using SDO data in the line A = 17.1 nm are displayed in 
Fig. 3(b). Apparently EF started moving before the flare, since its velocity was 
> lOOfcm/s one minute after the flare initiation. One can see that velocity of 
the earlier loop increased with time and exceeded that of the later loop; the 
first loop came up with the second one, thus leading to formation of the CME 
frontal structure. Note that beginning of the outer loop's motion was measured 
after beginning of the inner loop's motion. Till this moment, the outer loop 
was observed as a structure showing a weak brightness (see Fig. I(a2, b2)); 
its velocity was hard to determine. In what follows, the moment when we first 
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Figure 4. The same as in Fig. 1 but for CME on 13 June 2010. Unlike CME of 07 June 2011, 
initiation of CME in this case is related to eruption of flux rope whose nature is different from 
that of the filament (prominence). Label "1" on al, bl, b2 is the flux rope (outlined) before 
initiation of the CME-related flare, registered in the line A = 19.3 nm. 



determined translational velocity of the leading edge of its frontal structure (the 
CME front) will be referred to as the initial moment of propagation of the CME 
front. 

Fig. 3(c) shows time dependences of the CME velocity V(t) obtained using 
joint data from AIA, PROBA2 and LASCO, and soft X-ray intensity ISXR(t) 
from the CME-related flare region. Time profiles a(t) = dV(t)/dt and 
[dlsxn/dt](t) are presented in Fig. 3(d). Comparison between V(t) and IsxR,(t) 
profiles reveals that the time of increase in the CME velocity up to maximum 
value is comparable to the time of soft X-ray intensity enhancement (from the 
moment of the flare initiation to the moment of Isxr maximum); time of V(t) 
maximum is close to that of Isxn(t) maximum. These results agree with those 



in (Zhang and Dere, 2006 Maricic et al, 2007). Fig. 3(d) proves that time of 



a(t) maximum is close to that of [dlsxRl 'dt](t) maximum. This result agrees 

with that in (|Maricic et al, 2007[). 

" Temmer et al., 2010| ), time profile of the CME 



In (Temmer et al, 2008 



acceleration was compared to time variation in hard X-ray intensity lHx(t)- 
There were no sufficient data on hard X-ray emission from the flare region for 
that event; therefore, a(t) acceleration was compared to [dlsxii/dt](t). Such a 
change was possible due to the Neupert effect (Neupert, 1968). According to this 
effect, time variation in Inx(t) or in intensity of microwave radiation during the 
impulsive phase is close to dlsxa/dt time profile. 

Noteworthy is the fact that, after having reached its maximum value, the 
ejection velocity decreases fast (w 20 minutes) by several hundreds of km/ s; it 
thereupon passes to slow time variation. This influences also the time profile 



SOLA: paper.tex; 13 August 2012; 0:48; p. 8 



STUDY OF CME PROPERTIES USING HIGH RESOLUTION DATA 



600 




o l 1 0-1 , , 1 0,0 

5 10 1 10 100 

t, min(after flare onset - 05:32:52) t, min(after flare onset - 05:32:52) 

8, i 4-, , , ,4,0x10* 




1 , 1 -2-| , , — , , — , 1-1,0x10^ 

5 10 1 10 100 

t, min(after flare onset - 05:32:52) t, min(after flare onset - 05:32:52) 



Figure 5. The same as in Fig. 3 but for CME on 13 June 2010. Solid line in (a - b) denotes 
time dependence of velocity and of acceleration of the outer CME loop; dashed line, velocity 
and acceleration of the eruptive flux rope; dash-and-dot line, velocity and acceleration of the 
inner loop, according to SDO/AIA data. 



of the CME acceleration, see Fig. 3(d). Fig. 4 and 5 present CME recorded 
on 13 June 2010. This event has already been examined using SDO data in 
( Patsourakos, Vourlidas, and Stenborg, 2010). In what follows, we will present 
some new results. There are several important differences between this event and 
the previous one. Before initiation of the flare related to CME of 13 June 2010 
and before the CME formation, we observe activation and a slow motion of a 
loop- like structure at the base of the future CME source region (Fig. 4(al)). This 
is rather a loop of hot coronal matter than a filament, since it was not observed in 
H-alpha and Hcl 1083 nm lines during the period preceding the CME initiation. 
In order to distinguish between such loops and eruptive filaments (prominences) 
formed by colder chromospheric matter, the term 'eruptive flux rope' was coined. 
It emerged that images of CME of 13 June 2010 in lines at 17.1 and 21.1 nm 
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Figure 6. Difference images of CME of 08 March 2011, when CME is already formed, ac- 
cording to SDO data in the lines at 17.1 nm (a) and 21.1 nm (b); according to PROBA2 data 
(c); according to LASCO C2 data (d). (e) is the time dependence of the CME front velocity, 
according to joint data from SDO/AIA, PROBA2/SWAP and SOHO/LASCO (continuous 
dependence of velocity on time is denoted by solid line), and the time profile of soft X-ray 
intensity ISXR(t) - dashed line; (f) are the CME acceleration profiles (solid line) and dISXR/d 
t - dashed line. Vertical dashes show the time dependence of hard X-ray emission. 




t, min(after flare onset) 



Figure 7. (a) is the time dependence of the CME angular size 2a for four CMEs. ; (b) 
illustrates determination of the CME angular size 2a. Vertex of the angle determining the 
CME size is placed in the centre of the solar disc. 



SOLA: paper.tex; 13 August 2012; 0:48; p. 10 



STUDY OF CME PROPERTIES USING HIGH RESOLUTION DATA 




5 10 15 


<■■'<■ .,.l,i.'...-...»r 

20 


,,iil,,li,Vm,i;ll,l„,ir,rh,fi« 


(c) 




13.06.2010 ] 

. - i ii 



1.3 
1.2 
1.1 
10 



05:34 



05:36 



05:38 



05:40 



05:42 



Figure 8. (a) is the time variation in ratio of the CME height djj to its width dyy (square - 
13 June 2010; circle - 08 March 2011; triangle - 27 April 2011; diamond - 07 June 2011); (b) 
demonstrates the CME height and width; (c) is the comparison between [dn / dw\(t) (solid 
squares) and time dependence of 'aspect ratio'. Scale on the vertical axis on the left in (c) 
represents 'aspect ratio'; scale on the right, dn/dw 



differ significantly, as well as those of CME of 07 June 2011. This is evident from 
brightness distribution in the region of the frontal structure, see Fig. 4(c4). 

Like in the previous event, several high-brightness loops are observed at the 
initial phase of the CME formation (e.g., Fig. 4 (a3, b3)). In this case, inner loops 
catch up with outer ones too, thus forming the CME frontal structure. Occur- 
rence of the frontal structure is accompanied with occurrence of an extended 
dark area ('cavity'). 

Fig. 5 (a, b) shows that, unlike in the previous event, velocity of the outer 
loop (or the CME frontal structure) reaches its maximum value within 4 min- 



utes in the SDO/AIA field of view (as has already been found in (Patsourakos, 



Vourlidas, and Stenborg, 2010)); subsequently, the CME velocity decreases fast 



(within 7-9 minutes), too. Let us note one more difference between this event 
and the previous one: during the CME formation and its subsequent propagation, 
velocity of the eruptive flux rope is lower than that of the CME front in the 
greater part of its path. 

Time profiles of velocity and acceleration of the CME front (as deduced from 
joint SDO/AIA, PROBA2/SWAP and SOHO/LASCO data) are presented in 
Fig. 5 (c, d). During this event, like during the previous one, instants of time 
when V(t) and Isxn(t) as well as (t) and dlsxn/dt reached their maximum 
values are close. There is a slight difference between duration of the periods 
when V(t) and Isxn(t) reach their maximum values, too. 
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No solid line V(t) in Fig. 5 (c) is due to no reliable data on V(t) between the 
last measured location of the CME front (SDO/AIA data) and the first measured 
location of the CME front (the LASCO field of view) . Estimated average velocity 
in this interval of distances exceeds the maximum velocity value in Fig. 5 (a, c). 
This implies that there must be a hardly explicable maximum in this interval in 
V(t) profile. This problem can be solved if we assume, for instance, that CME 
propagates at a relatively large angle (« 30° — 45°) to the plane of the sky at 
the initial stage of its propagation, and that the angle between the CME axis 
and the plane of the sky decreases. Consequently, comparison between the CME 
velocities (i.e., projection of 3-D velocity onto the plane of the sky) measured in 
the AIA and LASCO fields of view becomes incorrect. 

The qualitative features of formation and propagation of CMEs of 07 June 
2011 and 13 June 2010 are typical of most events under consideration. The 
connection of their kinematic characteristics with soft X-ray intensity Isxn(t) 
from the flare region and with dIsxR./dt is typical, too. For the event of 27 March 
2011, we did not manage to establish connection between V(t) and Isxn(t), 
and a(t) and dlsxa/dt because of the complex Isxn(t) profile. Table 1 presents 
the main parameters describing kinematics of CMEs under consideration and 
its connection with soft X-ray emission from the flare region. For CMEs of 08 
March 2011 and 27 April 2011, decrease in velocity of the CME front V(t) after 
reaching its maximum value is relatively small (i.e., it is less than the arbitrary 
selected boundary of 100km /s). Let us exemplify it by the event of 08 March 
2011, Fig. 6 (e). Velocity V(t) decreases by ~ 95km/ s within 80 minutes, whereas 
the CME velocity on 07 June 2005, for instance, decreases by ss 570km/ s within 
37 minutes (Fig. 3(c)). Such variations in velocity with time influences the time 
profile of acceleration, see Fig. 6 (b): amplitude of absolute value of negative 
acceleration in this case is significantly lower than that of CMEs recorded on 13 
June 2010 and 07 June 2011. 



In( 


Fainshtein and Egorov, 2011 


Fainshtein and Zagainova, 2011 


). limb CMEs 


Fainshtein and Egorov, 2011 


) and halo CMEs ( 


Fainshtein and Zagainova, 2011 ) 



with insignificant variation of velocity V(t) after reaching the maximum value 
V max at moment t(V max ) form a special CME type. CMEs characterised by 
abrupt and dramatic variation V(t) at t > t(V max ) form another CME type in 
these papers. CMEs, whose velocity after reaching the maximum value varies 



insignificantly, agree with the CME propagation model proposed in ( Zhang and 



Dere, 2006). According to findings in (Temmer et al., 2008 Temmer et al. 



2010 Patsourakos, Vourlidas, and Stenborg, 2010 Fainshtein and Egorov, 2011 



Fainshtein and Zagainova, 2011 ) and in this paper, there are many CMEs whose 
velocity decreases fast by more than 100 km/s after having reached the maxi- 
mum. The CME kinematics model proposed in ( |Zhang and Dere, 2006 ) is not 



valid here. A more accurate model in this case is the model having 4 stages: 
there is also a stage of fast decrease in V(t) between the stages of the main 
acceleration and the insignificant variation of velocity. 

Correspondence between initiation of the CME-related flare and initiation of 
motion of the eruptive filament (flux rope having another nature) in the region of 
the CME initiation is of importance when selecting a possible mechanism of the 
CME generation. Correspondence between the flare initiation and beginning of 
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the main CME acceleration is of importance, too. Corresponding time is listed in 
Table 1. We can draw the following conclusions: on 27 April 2011, motion of the 
eruptive filament or flux rope precedes the flare; on 13 June 2010, 11 February 
2011, 27 March 2011 and 07 June 2011, it probably starts earlier. In these cases, 
determination accuracy of the initial moment of filament (flux rope) motion is 
insufficient. For the event of 08 March 2011, no eruptive filament or eruptive 
flux rope has been determined. As for the event of 27 April 2011, propagation 
of the CME front precedes the flare. It is possible that the standard flare model 
- CME (Temmer et al., 2010) - is valid for this event; according to this model, 
the CME propagation promotes initiation of a flare. 

Let us now discuss some general results concerning the initial stage of the 
CME propagation. 

1) For all events under consideration, the average duration of the main CME 
acceleration t acc 18.17 minutes was close to the average duration of enhance- 
ment of the soft X-ray intensity from the CME-associated flare region tsxR 
17.74 minutes (the event of 27 March 2011 was not taken into consideration in 
the second averaging). Here, t acc is the difference between the moment when 
CME velocity reaches its maximum V max and the moment when the first value 
of the CME velocity is determined from SDO data; t$xR is the duration of X-ray 
intensity enhancement Isxn(t), from the flare initiation to the maximum Isxr 
value. Besides, average values of the shortest duration of the main acceleration, 
combined with the shortest duration of soft X-ray enhancement (16.7 minutes), 
correspond to average values of the longest acceleration V max /t acc , combined 
with Vmax/tsxR. (0.98 km/s 2 ). Largest average values t acc combined with tsxR 
(27.6 minutes) correspond, conversely, to smallest average values V ma x/tacc com- 
bined with Vmax/tsxR (0.56 km/s 2 ). Both these results agree with results in 
( |Zhang and Dere, 2006[ ). 

2) . Dependences of the CME angular size 2a on time, as deduced from SDO 
and PROBA2 data, were plotted for 4 events (13 June 2010, 27 March 2011, 27 
April 2011, 07 June 2011), see Fig. 7(a). Fig. 7(b) shows the CME angular size 
as an angle with vertex at the centre of the solar disc and with rays traced to the 
remotest points (in latitude) of the CME boundary. For all CMEs in Fig. 7(a), 
2a increases with time. Range of variations in 2a is 1.5-3 for several ejections. 
The typical time of variation in the CME angular size (time of twofold variation 
in 2a since the first registration of CME) is w 3.5-11 minutes. 

3) . Peculiarities of time variations in ratio of the CME height dn to its width 
dw have been revealed for the same 4 events (see item 2): [djj /dw](t), Fig. 8(a). 
The CME height was counted from the CME source (arbitrarily considered as 
location of the CME-related flare, Fig. 8(b). During all the events presented in 
Fig. 8(a), dn/dw decreased with time either since the first moment for which we 
managed to calculate ratio dn / dw or since a later moment (see dependence for 
the event of 13 June 2010). This implies that the CME broadening exceeds its 
extension in the longitudinal direction during this period. For the three events 
(except that of 07 June 2011), a tendency to an insignificant variation in da /dw 
with time is observed. This can be interpreted as the transfer to the self-similar 
expansion of CMEs when shape of the CME boundary is not changed with time. 
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4. Discussion and conclusions 



Modern models of the CME formation are provided by a limited amount of 
experimental data on CMEs at early stages of their initiation and propagation 



(Forbes et al., 2006). The aim of this paper is to establish new regularities of 
the CME formation and of the initial stage of its propagation in order to use 
them for creating an adequate model of the CME initiation. This is possible due 
to data from instruments AIA aboard SDO that have unique characteristics, as 
well as due to joint data from PROBA2/SWAP, SOHO/LASCO and GOES. 

We have established physical and morphological features of formation and 
initial stage of propagation of six CMEs: 

1. event starts from the eruption of filament or loop-like structure having another 
nature. This has been reliably established for five events of 13 June 2010, 11 
February 2011, 27 March 2011, 27 April 2011, and 07 June 2011; 

2. soon after the filament eruption several overlying loops occur and start moving 
at different velocities, thus forming the CME frontal structure (for all the 
events under consideration); 

3. during the four events (13 June 2010, 11 February 2011, 27 March 2011, 
07 June 2011), velocity of the CME front decreases abruptly by more than 
100 km/s after having reached the maximum; it then varies insignificantly in 
the LASCO field of view. We suppose that such events form a special CME 
type. Velocity of two CMEs (08 March 2011 and 27 April 2011) decreases 
slightly after reaching the maximum value. Probably these CMEs represent 



another type of CMEs. According to the conclusion in ( Zhang and Dere, 2006 ), 
this pattern of velocity variation is typical for CMEs. Findings in (Temmer 
et al, 2008 Temmer et al, 2010 Patsourakos, Vourlidas, and Stenborg, 2010; 



Fainshtein and Egorov, 2011 Fainshtein and Zagainova, 2011) and in this 



paper prove that the most typical is the profile of the CME velocity where 
velocity decreases abruptly by more than 100 km/s after having reached the 
maximum. 

The formed CME frontal structure is not the same in different spectral lines. 
The line A = 17.1 nm is formed at a lower temperature than the line A = 
21.1 nm. The shift of boundaries of the CME frontal structure along the radius 
implies that temperature of outer regions of the CME frontal structure is higher 
than that of inner regions. 

Variation in the CME angular size 2a(t) demonstrates distribution and time 
variations in forces that influence the ejection boundary both outside and in- 
side CME. The conclusion concerning the increase in the CME angular size 2a 
between CME locations in the lower corona and in the LASCO field of view 
has been drawn in citePlunkctt. However, no data on typical temporal and 
spatial scales of such variations have been obtained. We managed to observe 
time variations in angular sizes of four CMEs in the SDO field of view. Value 2a 
was found to increase up to threefold at the initial stage of propagation. Ratio of 
the CME angular size in the LASCO field of view to the first measured value 2a 
could be 5 (for CME of 13 June 2010) and more. The typical temporal scale of 
doubling of the CME angular size after the first measurement of 2a was ss 3.5-11 
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minutes. In the frame of a simple CME model at the initial stage ('bubble', see 
( Patsourakos, Vourlidas, and Stenborg, 2010)), increase in the CME angular size 
with time evidences excess of internal force's influence on the CME boundary 
over external force's influence. 

Study of ratio of the limb CME's longitudinal size to its transverse size dn / dw 
provides solutions to several problems. First of all, this ratio is equivalent to 
the ratio of the CME velocity along its axis to the CME expansion velocity in 
the transverse direction. Knowledge of formation laws of this ratio allows us 
to estimate, for instance, propagation velocity of halo CMEs along their axes 
in three-dimensional space, by measuring the CME expansion velocity in the 
transverse direction (see discussion on this issue, based on LASCO data, in 
(Michalek, Gopalswamy, and Yashiro, 2009)). Besides, this ratio can be used 
to determine whether or not the CME expansion is self-similar. Comparison 
between CME propagation parameters and theory of self-similar expansion of 
CMEs ( Uralov, Grechnev, and Hudson, 2005 ) provides a more accurate estimate 
of self-similar CME propagation. If ratio du/dw remains constant over time, 
this implies conservation of the CME form (self-similarity) over time. 

Note that the parameter we have chosen to determine the CME expansion 
differs from that used in (Patsourakos, Vourlidas, and Stenborg, 2010 ) to perform 
the same analysis. That parameter was referred to as 'aspect ratio' (ratio of the 
height of centre of circle Hq , which outlined well most of the CME boundary, to 
the radius of this circle -Rf). Height He was measured relative to the solar limb. 



We used the parameter closest to that in (Michalek, Gopalswamy, and Yashiro, 



2009 ) . It can be shown that two these parameters should be connected by ratio 



Hq/Rf ~ 2g?h jdyif—X for the event of 13 June 2010. When deriving this formula, 
we ignored deviation of location of the CME-related flare from the limb and from 
the centre of the CME basis. Fig. 8 (c) compares dependence of He /Rf on time 
(Patsourakos, Vourlidas, and Stenborg, 2010) with our dependence [dn/dw] (t) 
for the event of 13 June 2010. We see that Hc/Rf variation over time is qual- 
itatively close to dependence [dn /dw](t)- In the first approximation, relation 
between the two parameters is described by their above-stated ratio. Noteworthy 
is the fact that, starting from instant t — 05:39:00, shape of the CME boundary 
docs not correspond precisely to circle (e.g., compare the visible CME boundary 



at t =05:40:11 in our figure with circle in Fig. 3 from (Patsourakos, Vourlidas, 



and Stenborg, 2010)). 



In conclusion, let us note that the use of data with high temporal and spatial 
resolution made it possible to obtain information about initiation and initial 
propagation stage of CMEs. However, a great deal needs to be done before 
initiation and development of CMEs become completely understood. In our next 
paper, we will try to relate the differences in V(t) profiles established for 2 CME 
types to the peculiarities of the magnetic field and other features of active regions 
wherein CMEs initiate. 
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